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We evaluate the decays B — > K*n, K(p,u, cf>) adding the long distance charming penguin contri- 
butions to the short distance: Tree+Penguin amplitudes. We estimate the imaginary part of the 
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parameterize its real part. The final results for branching ratios depend on only two real parameters 
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I. INTRODUCTION 



The CLEO II Q, Q , the BaBar 0, Q] and the Belle collaborations have recently reported data on the 
B meson non leptonic decay channels into a light vector and a pseudoscalar meson: 

B^PV. (1) 

These data are of the foremost importance for the determination of the angles of the unitarity triangle. 
In particular, if one of the two mesons in the final state is a strange particle the decays offer several 
channels for the extraction of the 7 angle, thus providing alternatives to the Ktt dec ay m ode. These non 
leptonic decays have been proposed long ago 0, B IE S Ell [HI Eil EH Edl 13 [ZU as a method for 
the extraction of the 7 angle. By this strategy one relates the data to a theoretical amplitude given as a 
sum of Tree and Penguin short distance contributions, whose relative weak phase is indeed 7. The usual 
computational scheme is the factorization hypothesis, either in the naive or in the QCD based version [l9| . 
20, 2lL l22u23l.l24| . Here one neglects non factorizable contributions as one proves that they are of the order 
of A/mb and therefore negligible in the mb — > 00 limit. 

There is a class of contributions however that, although suppressed in the infinite heavy quark mass limit, 
cannot be neglected a priori. These are long distance contributions known in the literature as charming 
penguin diagrams [25I I2H Eil EE Eof . Induced by the non leptonic hamiltonian oc V^V^bj^l — 
75)0 57^(1 — 75)s + Fierz + h.c, they cannot contribute in the vacuum saturation approximation. However, 
the insertion of at least two charmed particles (charm+anticharm) between the currents can produce sizeable 
contributions. As a matter of fact the 0(A x m^ 1 ) suppression is compensated by the Cabibbo-Kobayashi- 
Maskawa enhancement and the actual role of these contributions can be established only as a result of some 
dynamical calculations. In and we estimated these contributions for the B J PP decay channels 
and found that for the Ktt case they are indeed relevant. Their imaginary part can be evaluated with some 
confidence using the heavy meson chiral effective theory corrected for the light meson hard momenta. On 
the other hand, the real part is less predictable; for example in 0] and [2^ the results strongly depend on 
the cutoff in the loop integrals. 

In the present paper we give an estimate of the charming penguins in the charmless decays Q for the 
cases K*n and K(p,uj,(f>). We do not consider the particles 77,77' in the final state because the pseudoscalar 
SU(3) singlet is most probably contaminated by the glue and the evaluation of these contributions cannot 
be reliably done within our theoretical scheme. Due to the above mentioned difficulties we compute by 
the effective theory only the imaginary part of the amplitude, Tcv- For the real part of the amplitude, 
IZcv, we assume a simple parametrization IZcv oc lev, fixing the two proportionality constants (one for 
each of the SU(3)/ multiplets comprising the light vector mesons) by a fit. A similar calculation has been 
recently attempted by 30]; these authors assume a more phenomenological approach and parametrize the 
amplitude in the more general way (two complex numbers), using SU(3) flavor symmetry to relate the 
different amplitudes. Therefore our calculation can provide a dynamical test of the model in [30j . 

The outline of the paper is as follows. In Section [H] we estimate the tree diagram contribution (or 
Short-Distance contribution), computed in the factorization approximation. In Section ITTT1 we estimate the 
absorptive part of the charming penguin contribution. It is evaluated by use of the effective lagrangian 
for light and heavy mesons based on the Heavy Quark Effective Theory (for a review see [31^. The main 
uncertainty of this approach is the ex trap olation to hard light meson momenta and in this context we use an 
estimate of form factor given by us in 28] . In Scction HVI we present numerical results for the branching ratios 
and the asymmetry. In our previous papers we estimated the real part of the charming penguin contributions 
by using a cut-off Cottingham formula; the results where strongly cut-off dependent and for this reason we 
give here a simple parametrization of the real part and we estimate the relevant two real parameters by 
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comparing with the data. The overall result for both of the branching ratios and the CP-asymmetries is 
rather satisfying and points to a significant role of the charming penguins in the B — > K*7r, K(p, to, <j>) decay 
modes. 



II. SHORT-DISTANCE NONLEPTONIC WEAK MATRIX ELEMENTS 



The effective Hamiltonian for non-leptonic B decays is given by 

v: b v us ( Cl o^ + C2 o u 2 ) + v: b v cs { Cl oi + c 2 o c 2 ) - v; b v ts [ J2 ^ + Cg o g 



Ucs -T2 
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(2) 



where Cj are the Wilson coefficients evaluated at the normalization scale \x = mb [32], l33l l34l l35l |36( and 

next-to-leading order QCD radiative corrections are included. 0\ and 2 are the usual tree-level operators, 
Oi (i = 3, 10) are the penguin operators and O g is the chromomagnetic gluon operator. The Cj in eq. © 
are as follows 0: c 2 = 1.105, ci = -0.228, c 3 = 0.013, c 4 = -0.029, c 5 = 0.009, c 6 = -0.033, c 7 /a = 
0.005, cs/a = 0.060, cg/a = —1.283, cio/a = 0.266. Moreover, for the current quark masses we use the 
values 



m b = A.&GeV m u = AMeV m d = 8MeV m s = 0.150GeV. 
We define the T-matrix element 

(2tt) 4 S 4 { Pb - p w - p K , ) x Mk*-* = (K*n\ T \B) , 



(3) 



(4) 



with a similar definition for the K(p, u>, <p) final state. We separate short-distance and long-distance contri- 
butions to the weak B — > K*n decay: 



M = M S v + Mcv , 



(5) 



and evaluate the short-distance part of the amplitude M.sv using the operators in J2J) in factorization 
approximation (with Msv{B — > /) = —{f\TL c g\B) ). The results are in Tab. [I] Here ai = 



Ci+l 



Ci-l 



(i=odd) and cii = Ci -\ — — - (i=even). Moreover, if V is a vector meson, pW is a pseudoscalar meson we use 
3 

the following definition for the matrix elements of weak currents: 

(P(p)\A„\0) = -ifp Pll , (V(e,p)\V,\0)=f v m v e;, (6) 

and 



<PV)|vgp(p)) = Pi(g 2 



Fo(q 2 



m. 



(7) 



< V{e,p')\V» - A»\P{p) > 



2V{q 2 



M P + M v 
Mp + M v 



e^etpap'p - i(M P + Mv)e*Mi (q 2 ) 



2M- 



(p + p'YMq 2 ) + ^* ■ q) ~VV [A 3 ( q 2 ) - A ( q 2 )] , (8) 



where 



. , 2 , M v - M P A . 2 , 



M v + M P 
2Mv 



Ai(q 2 ) 



(9) 



and A 3 (0) = A Q (0). 
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TABLE I: Factorized Mst> amplitudes. 



Process 



Mst> 



B+ -> K* ^ 



+G F V2F^(m 2 K , ) f K . m K . V t * b V ts 



0,4 



aio 



B + -> K*+7T° 



-G F m K * {^f^Cm^O/if. [v* b V us a 2 - V t * b V t s(a 4 + aio) 



V^Vusai + V^Vts -(07 - 09) 



}(e*-Ps) 



B" - AT* *- 



V^Vus ai - Vj^Vts -(ag - 07) 



+F 1 B ^(m 2 K ,)f K . V* b V ts 



0,4 



O-10 



}(e*-PB) 



B° -» K* + 7T- 



- Gj? V2F^(m 2 K , ) /k. mK* KT 6 K S a 2 - V/ 6 Vt s (a 4 + aio) 



(s*-Pb) 



B+ -> A" V 



+ G F ^2A^ p {ra\) f K m„ V* b V ts 



(2a 6 - ag) m 2 K 



a 4 ~ — aio — 

2 (m b + m d )(m d + m s ) 



(s*-PB) 



B+ -» 



- G F rrip^Ay ~" p (m 2 K ) f K 
+ F 1 B ^(m 2 p )f p 



V* h V U s a,2 — VflVts 04 + aio — 



2(a 6 + as) m 2 K 
(m b + m u )(m u + m s ) 

V^Vu, 01 - I VJ* fc Vi. (07 + 09)] } (e* • ps) 



- G F m p ^F^ K \m 2 p ) f p V* b Vus ai - K ( * 6 V ts ^(a 7 + 09)] 

(2a 6 - as) rn 



+ A^ p (m 2 K )f K V t * b Vts 



2 

a4 — ^ ai ° ~~ 7 w ~ — 7 

2 (m b + m d )(m d + m s ) 



} (e* 'Ps) 



-G i ,V2A 3 ^(m| f )/ Jf m, 



/ 2(a6 + as) m?, \ 



(e*-PB) 



B+ -» A"+w 



-Gf {^-^(m*)/,., 



V* b V us 02 — V t * b Vts [ 0,4 + am 



- i(a 7 +ag)j 




2(a6 + as) 





(m b + m u )(m u + m s ) 



(e* ' Ps) 



-Gjr m w jif^ (m* ) / u l^V^oi - V t ;Vt s ^2(03 + a 5 ) + i(a 7 + 09)) 

(2a6 - as)m ; 



- A B ^(m 2 K ) f K V; b V t s a 4 - -aio 



(m b + m d )(m d + m s ) 



> (e* -Pb) 



+G F ^/2F^ K {m 2 ) U V t * b V t 



a 3 + a 4 + 05 — 



07 + ag + aio 



(e*-Pfl) 
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III. ABSORPTIVE PART 



The computation of the discontinuity of the charming penguins diagrams contributing to B — > K*ir gives 
(cfr. analogous diagrams in figure 2 of Ref. 

j4. 



/d 
(27Tj 4 



S+(q ~ m D )S+(p DM - m D )x 



x M(B-> dWdW)M(dMdM -> K*n) = 



y/uj* 2 - 1 J dnM(B -» DWflW) M(I)WdW ^ X*tt) 



... v ~ - , . ~ ... • (10) 

where the integration is over the solid angle and a sum over polarizations is understood. A similar equation 
holds for the K(p,u>,(f>) final state. The amplitudes for the decays B — ► DgD^ are computed by factor- 
ization and using information on the Isgur-Wise function (see below). Diagrams for the calculation of the 
amplitudes — > K*n are in FigJI] A similar diagram can be drawn for the Ds ~ * K(p,u>,4>) 
amplitudes. 



K 



7t 



K 



K 



(a) 



(b) 



(c) 



K 



71 



K 



71 



(d) 



(e) 



FIG. 1: Diagrams for the calculation of the amplitudes 7J'*' — > K*ir amplitudes 



The effective lagrangian to compute these diagrams can be written as follows (see e.g. |31| ) 

L = i < HbifDfojEta > +W_< Hb^jsA^Ha > 

+ i(3 < H b v» (V„ - p^) ba H a > +zA < H b a^F^(p) ba H a > . 



Here < ... > means the trace, 



(11) 



(12) 



and P5, P M are annihilation operators of the pseudoscalar P and vector V heavy mesons normalized as 
follows: 



\Pk\P)=VMH, 



\P„\V(e)) = e^y/M^ 



(13) 
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The first term in the lagrangian (|llfl contains the kinetic term for the heavy mesons giving the P and P* 
propagators, i/(2v ■ k) and — ?(<? Ml ' — v tJ, v u )/(2v ■ k) respectively. However, since the residual momentum is 
not small we will use the complete form of the propagators instead of their approximate expressions [3, [2^] . 
The interactions among heavy and light mesons are obtained by the other three terms. In the second term 
there are interactions among the heavy mesons and an odd number of pions coming from the expansion A^ 
. Here A^ — \ (<^<9 M £ — ^6^). Expanding the field £ = exp(z7r//) and taking the traces it provides the 
coupling in the lower vertices of Fig. ^ The last two terms give the upper vertices of Fig. ^ The octet of 
vector resonances (p. K* , etc.) is introduced as the gauge multiplet according to the hidden gauge symmetry 
approach of Ref . [37J • We put 

where p is the usual hermitian 3x3 matrix of flavor SU (3) compri sing the nonet of light vector mesons and 
gv is determined by vector meson dominance as follows: gv — 5.8 |37j ]. 

Also the parameter can be fixed by vector meson dominance. This corresponds to assume that, for the 
heavy pseudoscalar mesons H , the coupling of the electromagnetic current to the light quarks is dominated 
by the p, u>, 4> vector mesons. This produces the numerical result 

/3 = ^«0.9. (15) 

gvjv 

For q we take the recent experimental result from the CLEO collaboration, obtained by the full width of 
D*+ 38]; this gives a value: g = 0.59 ± 0.07 ± 0.01. 

It can be noted that while these determinations consider soft pions we are interested in the coupling of a 
hard pion to D and D* . This introduces a correction that can be parametrized by a form factor, see below 
for a discussion. 



Let us now consider the parameter A in Eq. (|11|) . In |31| one can found an estimate of this parameter 
based on the effective chiral lagrangian for heavy mesons. We present here an update of this analysis and 
new numerical results. 

It is useful to begin by writing down the phenomenological heavy-to- heavy current in the leading 1/toq 
approximation: 

Jt = ■ v ') < ^ c) 7m(1 ~ 75)^ fc) > , (16) 

which gives rise to 



< D(v')\c^b\B(v) > = y/M B M D £(« • v') (v + «% , 
<D*{v',e)\c lll b\B{v)> = y/M B M D £(v ■ v') e^eVt/" 



< D*(v', e)\c ltil5 b\B{v) > = y/M B M D £(v ■ v') i[(l + v i/)ej - (e* • v)t£] , (17) 

where £(v ■ v') is the Isgur-Wise form factor. Note that, by this definition, the choice of phases in Eq. i|17|) 
agrees with © and 10. We also write the phenomenological heavy-to- light leading current for coupling to 
light pseudoscalar mesons of the heavy mesons in the multiplet H 

K= Y<7 M (1- 75)^4 a > ■ (18) 
Here F is related to leptonic decay constant of the pseudoscalar heavy meson B appearing in © by 



(19) 
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Several numerical analyses based on QCD sum rules have been performed, see for a review e.g. |3l|. The 
result we shall use is 

F = 0.30 ± 0.05 GeV 3/2 (20) 

which is obtained neglecting radiative corrections since they are neglected also in the evaluation of the 
Isgur-Wise form factor the parameterization we use below is based on. 

The effective theory approach gives predictions for the form factor V(q 2 ) at high q 2 and relates it to 
the parameter A. We can match this result with the theoretical calculations coming from Light Cone 
Sum Rules 0, (LCSR) and lattice QCD 0, 0- To do this we compute the form factor at 
q 2 ~ q^nax where it is dominated by the nearest low-lying vector meson pole. Computing the form factor at 
q 2 ~ q 2 max = {Mb — My ) 2 and at leading order in I/tuq, one gets 0: 

vi 2 ^ gv^ Mp + M v 1 

V(q^) = ^XF v _ — - (21) 

where A is the appropriate mass splitting, i.e., for the transition B — > K*, A = tub* —vtib — 135 MeV. From 
the LCSR and lattice QCD analyses for the transition B — > K* we infer the value V = 1.5 at q 2 = q 2 = 17 
GeV 2 . In order to compare with Eq. I|21|l we assume that in the range q 2 S (q 2 , q^^) the form factors are 
dominated by the B* simple pole. This gives V{q 2 nax ) = 1.8 and, as a result, 

A = +0.56 GeV" 1 . (22) 

We note that to determine A we used a positive sign of V(q 2 nax ). The same sign would be obtained using 
HQET and by assuming that the strange mass as large, see eq. I|17|l . We also note that a different result for 
the phase and the magnitude of A was obtained in |3lJ . The difference in magnitude was due to the different 
methods used. Here we use LCSR and lattice QCD, while in [il| an extrapolation is used from the D — > K* 
form factor. As to the phase, it is due to a different choice we use here for the phase factor of the heavy 
vector meson, see eq. (|12ll . 

To compute the matrix elements we use the following kinematics: 

p^ = m B i/' ■= (m s ,0) , p£ ( „, = m D v'^ , q = p-p DM . (23) 
m 2 B + m 2 D — m 2 D 

where to* — and the angular integration is over the directions of the vector v = 

ImB-rriB 

nVu>* 2 — 1. Our results are as follows: 

M(B(v) -» D s (q)D*(e,v')) = -K (m B + m D ) e* • v , (24) 
M(B(v) — > D*(i],q)D(v')) = -K m Ds r)* ■ (v + v') , (25) 
M(B(v)^D* s (ri,q)D*(e,v')) = - iKm D ^t* a (ie aX ^v ,x v c - g m (l + uj*) + v a v'^) , (26) 

G F 

where K = —=V* b V cs a2y/mBTnDfD s £iw(w*) ■ On the other hand, for the K*ir final state we have 
v2 

M(D s {q)D {e,v')->K {p K , e)Tr(p n )) = -=J e x e a 

Jrr V2 V m D s 



2l3m D q°pl AXm D ,G aX {p^,p K ,v') 



(m D v ; - Pk) 2 - m 2 D (m D *v' - p^) 2 - m\ 



(27) 
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M(D* s ( V ,q)D(v') -> K*(p K ,e)n(p^)) = - 



2gm D , F 2 (\pk\) gv I m D 



V% V m Ds (m D *v' -Pk) 2 



D' 



2pq a (+pl - ^-^p\] -A\mn a H° X (j>*,PK,v') 



(28) 



Jtt v 2 V 



4Am J >g a (p J f) A1 pg 5181 

(m D v' - Pk) 2 - m 2 D m D * (m D *v' - p^) 2 - m 2 D , 



{2[3q°5 T x + A\m D ,{p T K 5l - ( PK )xg aT )) ,(29) 



where 



G a \p^p K ,v')= - W ■q)^g aX (pK-p,)-p:p X K^j - (q-p,)(v'°p x K -g a \v' -pk)^ 
- q X (p°(pK-v')-v' a ( PK -p^ , 



(30) 



H° x (p„,p K ,v') = g^(p K .p^-^( m ^- PK . q ))-p^ K (p^ + ^LqA . (31) 

Eq. I|27|> corresponds to the sum of diagrams (a) and (b) of Fig. ^ Eq. (|28|) corresponds to diagram (c) 
and, finally, Eq. (|29|) corresponds to the sum of diagrams (d) and (e). Similar results hold for the K(p,u, <fi) 
final states and we do not reproduce them here for brevity sake, -F(|p,r|) is a form factor taking into account 
that in the vertex DD*tt the pion is not soft and therefore the coupling constant should be corrected. Its 
determination by a quark potential model is discussed in |2^| using a quark model. The central value we use 
is FdpVI) = 0.065. In absence of more detailed information, the same form factor is adopted for the upper 
vertices in Fig. as well. 



IV. NUMERICAL RESULTS 



Let us first consider the short distance contribution. Using, for the CKM matrix elements, p — 0.2296, 
r\ = 0.3249, A = 0.819, and for the involved form factors the value collected in Table IV of Ref. |45j|. one 
gets the results reported in Tab. |H] These values correspond to the following value of the angle 7 of the 
unitarity triangle: 

7 = arctan f ^ ~ 54.8° . (32) 

Next we consider the long distance absorptive part. For its parameters we assume the values of previous 
Section. For the D s decay constant we use fn, = F / Jmu B , with F given in l|20l) ; for the Isgur-Wise form 
factor we use the parameterization 

ec^-^) = f-— ^ — 7 V (33) 

\1 + v ■ v' J 

with p 2 ~ 1. This parametrization agrees with the results for the b —> c exclusive decays obtained by the 
CLEO and Belle Collaboration. It is a fit to the results obtained by the QCD sum rules method, see |3J. 
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TABLE II: Theoretical values for Mst>, Mct>- Units are GeV. 



Process 


A4 5 -d x 10 s 


lmMcv x 10 8 


B+ -> K*°n+ 


+ 1.45 


-2.14 


B+ K*+n° 


+ 1.02 - 0.79 « 


-1.52 


B° -» A'*%° 


- 0.60 - 0.08 i 


+ 1.52 


B° -» K'+n- 


+ 0.85 - 1.01 i 


-2.14 


B + _> A"V 


+0.17 


-2.74 




+ 0.39 - 0.64i 


-1.94 


B° -» A>° 


+ 0.48- 0.11 i 


+1.94 




- 0.28 - 0.74i 


-2.74 


B+ -> AT+o; 


-0.27 - 0.63 i 


-1.94 


B° -> #°<j 


+ 0.03-0.10i 


-1.94 


B+ -» A+0 


1.30 


-2.58 


B° -» A"V 


1.30 


-2.58 



The numerical results for the imaginary part of M.cv are given in Tab. HH Typical sizes of the different 
contributions to B~ — > A' 7r~ arc as follows. From the two terms in Eq.J37|): lmA4%j = —1.32 x 10 -8 ; from 
the terms in Eq.|2EJ|: lmM c cv = -0.16 x 10~ 8 ; from the terms in Ea.QSfy: hxiM.% = -0.66 x 10~ 8 . For the 
B~ -> K°p- channel wc find: ImM^ = -1.52 x 10~ 8 ; hnM c cv = -0.70 x 10" 8 ; ImM% = -0.52 xl0~ 8 . 
It can be noted that the phase of M.sv is purely weak while the phase in Mcv is only due to strong 
interactions. 



A. Branching ratios 

From the results of Tab. [H] we can compute the Branching Ratios (£>) and the CP asymmetries. As 
explained in the Introduction, we have not presented an evaluation of the real part of Acv- Attempts for 
the K ix channel can be found in [3, . Typical results are that the real part is of the same order of the 
imaginary part, but the uncertainties are large. To estimate the CP averaged Branching Ratios we add 
therefore to the imaginary part a real part as follows: 

M cv — TZcv + ilcv (34) 

and consider the values TZcv — 0, ±Tcv- The results are reported in Tab. IIIII for the three values of IZcv, 
together with the results obtained considering only the short distance amplitude. They show the relevance 
of the long-distance contribution. The "best value" of TZcv might be obtained by a fit to the available 
experimental data; these data are reported in Tab. IIVI and a comparison with the results of the previous 
table shows that the preferred values, except for the channels with <fr in the final state, satisfy TZcv ~ ~~^cv- 
Instead, the decay in K 4> prefer TZcv ~ +%cv- To find a "best value" of TZcv we defined a \ 2 as: 

where the Br^ xp are reported in Tab. IIVI and the a(Bri) are obtained summing quadratically the statistic 
and systematic errors in the same table. When the experiment provides an asymmetric error a conservative 
error was assumed: a = Max(a\, 02). For the decay B° — > K*°ir° we put Br exp = and the corresponding 
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TABLE III: Branching Ratios B (units 10 -6 ). In the second and third columns, theoretical values computed using 
only the short distance amplitude and the full amplitude (i.e. short distance and long distance); in the latter case 
the three values correspond respectively to 1Zct> = (+1, 0, —1) x Xct>, see text. 



Process 


B (SV only) 


B (SV+CV) 


B + K*%+ 


1.96 


(4.71,6.22,16.3) 


B + _> K*+7v° 


1.56 


(5.20,5.95,11.0) 


B° -> K*°n° 


0.32 


(2.51,2.10,5.66) 


B° -» K*+tt- 


1.50 


(9.94,9.13,16.2) 


B + _» K°p+ 


0.03 


(13.1,6.99, 14.8) 


B+ _» K+p° 


0.52 


(8.42,6.32,11.2) 


B° -> A°p° 


0.21 


(7.88,3.07,4.71) 


B° - A+p" 


0.54 


(18.1,10.4,15.6) 




0.44 


(10.7,6.21,8.73) 


B° -» A°cj 


0.01 


(6.70,3.58,6.91) 


B+ -> 


1.55 


(7.58, 7.63,19.8) 


B° -> K°<f> 


1.42 


(6.98, 7.02,18.3) 



error was fixed to be equal to the experimental upper limit a = 3.6 x 10 . We have first attempted a fit 
with only one real parameter r = TZcv /lev , which corresponds to the use of the SU(3) nonet symmetry for 
the light vector mesons. The minimum value of x 2 (= 12.1) is obtained for r = 1.207. Since however we do 
not expect the validity of the nonet hypothesis we have also tried a two parameter fit, corresponding to the 
two multiplcts (octct+singlct) of the light vector mesons. This gives as a result TZcv = — 0-954Tcv for the 
(p, K*, ui) set of particles and two solutions, TZcv = — 0.2012£d and TZcv = +l-21Xcv, for the <f> with the 
same \ 2 — 7-8. The former solution would be preferred because it represents a less important deviation of 
the SU(3) nonet symmetry. Even if the two solutions produce the same values of the Branching Ratios they 
would produce different CP asymmetries. The results we find are in Tab. IIVI Besides our data, we also 
present two model calculations presented in [30j . where an analysis of all the PV (not only strange) final 
states is performed. The first model (called in that paper Scenario 1) contains short distance terms (QCD 
factorization) and unconstrained 7 angle. Other theoretical determinations based on QCD factorization are 
in Ref. Ilcl. e.g. B(B° — * K + p~) = 12.1. The authors of jstj do not agree with the conclusions of Du et 
al., Ref. [46j, because, differently from them, they include the K*tt channels. The second model (Scenario 
2) uses QCD factorization and charming penguins as in the present paper, but differently from this paper, 
where only two real free parameters are used, they employ two universal complex amplitudes multiplied by 
a computed Clebsh-Gordan coefficient. It is worthwhile to note that a significant agreement with the data 
is obtained with our simple hypothesis. We also note that the value 7 of Eq. (|32f> obtained by a global fit 
to the CKM matrix 01 is compatible with the B — > K*ir ,K(p, u>, 4>) Branching Ratios only if the charming 
penguin diagrams are included. 
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TABLE IV: CP averaged Branching Ratios B (units 10 -6 ). In the second column theoretical values computed using 
the present model. In third column and fourth column theoretical computations based on Ref. |3o| : Scenario (Sc.) 1 
refers to QCD with factorization and free 7; Scenario 2 refers to QCD+Charming penguins with constrained 7 (see 
text). Experimental data are from CLEO 0,0, BaBar 0,0, and Belle or averages from these data. 



Process 


6 (this paper) 


B ([22], Sc. 1) 


B ([22], Sc. 2) 


B(Exp.) 


B+ -» A'*°7T+ 


15.6 


7.889 


11.080 


12.1 ±3.1 (av.) 


B~ -> A*-7T° 


8.44 


7.303 


8.292 


7.r^ 1 1 4 ± 1.0 (< 31) (CLEO) 


B° -» K*°n° 


5.61 






< 3.6 (CLEO) 


B° -» K*+7v~ 


12.0 


9.760 


10.787 


19.3 ±5.2 (av.) 


B + -> K°p+ 


14.1 


7.140 


14.006 




B + K+p° 


8.56 


1.882 


5.665 


8.9 ±3.6 (av.) 


B A°p° 


4.86 


5.865 


8.893 




S ° _> A+p" 


11.6 


6.531 


14.304 


15.9 ±4.7 (av.) 


B+ -» A+lj 


6.19 


2.398 


6.320 


9.2l|j ± 1.0 (CLEO); < 4 (BaBar); < 7.9 (Belle) 


B° -» A'°lj 


6.27 


2.318 


5.606 


6.3 ±1.8 (av.) 


B+ -> A+0 


9.11 


8.941 


9.479 


8.9 ± 1.2 (av.) 


B° -> AV 


8.39 


8.360 


8.898 


8.7 ±1.4 (av.) 



B. CP asymmetries 



From previous results we can also compute the CP asymmetries for the various channels. The Belle 
Collaboration reported the result 



A K 



B(B~ -v K~u>) - B(B+ -► K+lu) 



CP 



B(B~ -> K~u) 
The result we find agrees with the data: 



K+L 



-0.21 ±0.28 ±0.03 . 



A 



CP 



= -0.37 



(theory) . 



We also get 



A 



go _ B(B° -> uK°) - B(B° -> ljK ) 
CP ~ B{B° cjK°) + B(B a -> ujK°) 



= -0.05. 



We can similarly compute the asymmetries for the AT*7r, Kp channels, defined by 

At 

■At 
A 00 



B(B+ - 


-> K*+ir°) 


-B(B- 


-> A-*-7T°) 


B(B+ - 


-4 A"* + 7T°) 


+ B{B- 


-> A^*-7T ) 


B(B+ - 


-> K*°ir+) 


- B(B~ 


— *0 . 
— > A 7T ) 


B(B+ - 


-> A'*°7T+) 


+ B(B- 


-> A 7T~) 


B(B°- 






-»■ A~*~7T+) 


B{B° - 


* A"* + 7r^) 




-> AT*-7r+) 


B{B° - 


+ A"*°7r°) - 




> ;TV) 


B(B° - 


+ A"*°7T ) 4 




> iTV) ' 



(36) 



(37) 



(38) 



(39) 



with analogous definitions for the Kp channel, with A^ — > *4 P and the changes A"* — ► if and 7r — > p. 
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TABLE V: Theoretical values for the asymmetries; for the definitions see 13911 . 



Asymmetry 


A +0 


A 0+ 


A + ~ 


.4 00 


.4+° 


.4°+ 


A + ~ 


A 00 


A C p 


0.27 





0.31 


-0.04 


0.27 





0.30 


-0.08 



The results are presented in Table They have a peculiar pattern and will therefore provide a crucial 
test of the present model when future experimental data for these asymmetries are available. 
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